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Dark Matter (DM) is a fundamental ingredient of our Universe and of structure formation, and 
yet its nature is elusive to astrophysical probes. Information on the nature and physical properties 
of the WIMP (neutralino) DM (the leading candidate for a cosmologically relevant DM) can be 
obtained by studying the astrophysical signals of their annihilation/decay. Among the various e.m. 
signals, secondary electrons produced by neutralino annihilation generate synchrotron emission 
in the magnetized atmosphere of galaxy clusters and galaxies which could be observed as a diffuse 
radio emission (halo or haze) centered on the DM halo. A deep search for DM radio emission 
with SKA in local dwarf galaxies, galaxy regions with low star formation and galaxy clusters (with 
offset DM-baryonic distribution, like e.g. the Bullet cluster) can be very effective in constraining 
the neutralino mass, composition and annihilation cross-section. For the case of a dwarf galaxy, 
like e.g. Draco, the constraints on the DM annihilation cross-section obtainable with SKAl-MID 
will be at least a factor ^10^ more stringent than the limits obtained by Fermi-LAT in the y- 
rays. These limits scale with the value of the B held, and the SKA will have the capability to 
determine simultaneously both the magnetic held in the DM-dominated structures and the DM 
particle properties. The optimal frequency band for detecting the DM-induced radio emission is 
around ^ 1 GHz, with the SKAl-MID Band 1 and 4 important to probe the synchrotron spectral 
curvature at low-v (sensitive to DM composition) and at high-v (sensitive to DM mass). 
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1. Unveiling the nature of Dark Matter 

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned, and 
it is a unique multi-disciplinary experiment. Even though the SKA, in its original conception, 
has been dedicated to constrain the fundamental physics aspects on dark energy, gravitation and 
magnetism, much more scientific investigation could be done with its configuration: the exploration 
of the nature of Dark Matter is one of the most important additional scientific fhemes. 

Among fhe viable compefifors for having a cosmologically relevanf DM species, fhe leading 
candidate is fhe lighfesf particle of fhe minimal supersymmefric exfension of fhe Sfandard Model 
(MSSM, Jungman ef al. 1996), plausibly fhe neufralino %, wifh a mass in fhe range befween 
a few GeV fo several TeV. Information on fhe nafure and physical properties of fhe neufralino DM 
can be obfained by sfudying fhe asfrophysical signals of fheir inleracfion/annihilalion in fhe halos 
of cosmic sfrucfures. These signals (see Colafrancesco ef al. 2006, 2007 for defails) involve, in fhe 
case of a DM, emission along a wide range of frequencies, from radio fo y-rays (see Fig. for 
a DM specfral energy disfribufion (SED) in a fypical dwarf galaxy). Neufral pions produced in XX 
annihilation decay prompfly in -A yy and generate mosf of fhe confinuum phofon specfrum al 
energies E>\ GeV. Secondary elecfrons are produced fhrough various prompf generafion mecha¬ 
nisms and by fhe decay of charged pions, ^ + v^(v^), wifh ^ + v^(v^) -|- Ve(Ve). 

The differenf composition of fhe XX annihilalion final slate will in general affecl fhe form of fhe 
elecfron specfrum. The time evolution of fhe secondary election specfrum is described by fhe 
fransporl equalion: 

^ = V [DVn,] + ^ [be{E)n,] + Qe{E,r) , (1.1) 

where Qg{E,r) is fhe source specfrum, ng{E,r) is fhe equilibrium specfrum (al each fixed 
lime) and bg = —dE/dt is fhe energy loss per unif lime, bg = bjcs + bsynch + btrem + bcoui (see 
Colafrancesco ef al. 2006 for defails). The diffusion coefficienl D in eq.(|T^ sels fhe amounf of 
spalial diffusion for fhe secondary elections: if lurns oul lhal diffusion can be neglecfed in galaxy 
clusters while if is relevanf on galactic and sub-galaclic scales (see discussion in Colafrancesco ef 
al. 2006, 2007). Under fhe assumption fhal fhe population of high-energy can be described by a 
quasi-sfalionary {dtigjdt « 0) fransporl equation, fhe secondary election specfrum ng{E,r) reaches 
ifs equilibrium configuration mainly due fo synchrolron and ICS losses al energies E > 150 MeV 
and due fo Coulomb losses al lower energies. Secondary elections evenlually produce radialion 
by synchrolron emission in fhe magnelized almosphere of cosmic sfrucfures, bremssfrahlung wifh 
ambienl protons and ions, and ICS of CMB (and olher background) photons (and hence an SZ 
effecl, Colafrancesco 2004). These secondary parlicles also heal fhe ambienl gas by Coulomb 
collisions. 

A large amounf of efforls have been pul in fhe search for DM indirecl signals al y-ray en¬ 
ergies looking predominanlly for Iwo key specfral feafures: fhe -A yy decay specfral bump , 
and fhe direcl XX 77 annihilalion line emission, wifh resulls lhaf are nol conclusive yel (e.g., 
Daylan el al. 2014, Weniger 2012, Doro el al. 2014). The non-defection of signals related to DM 
annihilation/decay from various asfrophysical largefs (including observalions of dwarf spheroidal 
galaxies, fhe Galactic Cenler, galaxy clusfers, fhe diffuse gamma-ray background emission) is in- 
lerpreled in terms of conslrainfs on fhe (self-)annihilalion cross-seclion (or decay lime) of fhe DM 
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Figure 1: Left. The Draco dSph multi-wavelength spectrum for a 100 GeV WIMP annihilating into bb. 
Right. The effect of varying the magnetic field strength on the Draco multi-wavelength spectrum for a 100 
GeV WIMP annihilating into bb. The WIMP pair annihilation rate has been tuned as to give a y-ray signal 
at the level of the EGRET measured flux upper limit (from Colafrancesco et al. 2007). 


particle candidate. Assuming, for instance, a canonical WIMP of = 100 GeV annihilating to 
Z^-quarks, stacked observations of dwarf spheroidal galaxies with Fermi-LAT put a constraint of 
(aV) < 2 X lO^^^cm^s^^ (95% c.l.) on the thermally-averaged self-annihilation cross-section of 
DM particles (Ackerman et al. 2014). Hopes of discovering annihilating WIMPs in /-rays are 
relegated to Fermi-LAT successors and the forthcoming CTA experiment (Doro et al. 2013). 

There are, however, also good hopes to obtain relevant information on the nature of DM from radio 
observations of DM halos on large scales, i.e. from dwarf galaxies to clusters of galaxies. 

2. Radio emission signals from Dark Matter annihilation 

Observations of radio halos produced by DM annihilations are, in principle, very effective 
in constraining the neutralino mass and composition (see, e.g., Colafrancesco and Mele 2001, 
Colafrancesco et al. 2006, 2007), under the hypothesis that DM annihilation provides an observable 
contribution to the radio-halo flux. 

The wide range of frequencies probed by the SKA and the variety of achievable observational 
targets (and in turn of magnetic fields) will allow testing the non-thermal electron spectrum from 
about 1 GeV to few hundreds of GeV, that is the most relevant range in the WIMP search. Fig¬ 
ure ^ displays the predicted spectral differences between various annihilation channels and WIMP 
masses: note that these differences manifest mainly in low-v slope variation for differing annihila¬ 
tion channels and in the high-v spectral flattening/steepening for larger/smaller masses. The SKA 
sensitivity curves for SKAl-LOW and SKAl-MID are taken from Dewdney et al. (2012). 

The surface brigthness produced by DM-induced synchrotron emission is heavily affected by dif- 
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Figure 2: Flux densities for dwarf galaxies (M = 10^ Mq), galaxies (M = 10^^ M©), and galaxy clusters 
(M — 10*^ Mq). The halo prohle is NFW, {B) — 5 /rG, and the annihilation cross section is hxed to the 
value (cry) ~ 3 x 10^^^ cm ^ s^^ Black lines are for dwarf galaxies, red lines for galaxies and green for 
clusters. Left; WIMP mass is 60 GeV, solid lines are for composition bb and dashed lines for T+T^. Right; 
Composition is bb, solid lines are for WIMP mass 60 GeV and dashed lines for 500 GeV. Dash-dotted lines 
are SKA sensitivity limits for integration times of 30, 240 and 1000 hours (Dewdney et al. 2012). The flux 
is calculated within the virial radius (from Colafrancesco et al. 2014). 

fusion in small scale structures, e.g., dwarf and standard galaxies, while it is less important in large 
structures, e.g., galaxy clusters (see Colafrancesco et al. 2006-2007 for a detailed discussion). 

Polarization from DM-induced radio emission is expected at very low fractional levels due to 
the fact that DM spatial and velocity distribution is nearly homogeneous and that DM annihilation 
is mediated by secondary particle production. Therefore, a low polarization level of detectable 
radio signals in the directions of DM halos would be consistent with the DM origin of such radio 
emission. Residual high-polarization signals could be hence attributed to astrophysical sources in 
the direction or within the DM halos, and one could use these signals to infer properties of the 
magnetic field in these structures (see R. Beck et al. 2015, and F. Govoni et al. 2015). 

2.1 Cosmological evolution of Dark Matter radio emission 

Figure ^ displays the evolution of radio emission from DM halos of mass 10^ Mq, 10^^ Mq, 
and lO'^ Mq for a constant magnetic field of 5 /fG, in accordance with arguments made in Co¬ 
lafrancesco et al. (2014). The annihilation channel is bb, the mass of the neutralino is 60 GeV and 
a DM annihilation cross-section {oV) = 3 x 10^^^ cm ^ s^* was adopted. Emission from dwarf 
galaxy halos {M k, 10^ M©) is just below the SKA detection thi'eshold for this value of {oV) but 
would be visible at redshift z < 0.01 for the assumed DM annihilation cross-section. This justify 
the search of DM-induced radio signals mainly in dwarf galaxies of the local environment, at dis¬ 
tances ^ 3 Mpc (z 0.0007). Emission from galactic DM halos (M R:; 10*^ M©) are detectable 
by SKA out to z ~ 0.8 even with the reference value of {oV), and can provide a non-detection 
upper-bound on {oV) an order of magnitude below the assumed value even at such high redshifts. 
Emission form galaxy cluster halos can provide similar constraints but out to higher redshifts z ^ 3. 
These objects thus offer the option of deep-held observations that can scan a larger fraction of the 
DM parameter space than the best current data. 
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Figure ^ also shows that the effect of diffusion are far less significant when observing higher-z 
objects, again simplifying the modelling and analysis of the SKA observations. 



Figure 3: Flux densities for various DM halos within the virial radius. Left; dwarf spheroidal galaxies 
(M = 10^ M 0 ); mid panel: galactic halos (M — 10^^ M©); right panel: galaxy clusters (M = lO'^ M©). 
We assume a NFW halo profile, (cjy) = 3 x 10^^^ cm ^ s^^ and {B) = 5 pG. WIMP mass is 60 GeV and 
the composition is bb. Solid lines are without diffusion, and dotted are with diffusion. The SKA sensitivity 
limits (dash-dotted lines, Dewdney et al. 2012) are shown for integration times of 30, 240 and 1000 hours, 
respectively (from Colafrancesco et al. 2014). 


2.2 Optimal DM laboratories 

In order to identify the optimal DM laboratories for radio observations we scan a parame¬ 
ter space extending from dwarf galaxies to galaxy clusters over a wide redshifts range z « 0 — 5. 
The choice to examine the halos of both large and small structures is crucial, as dwarf spheroidal 
galaxies are well known to be highly DM dominated but produce faint emissions, while larger struc¬ 
tures, but not immaculate test-beds for DM emissions, provide substantially stronger fluxes. This 
indicates that a survey of DM halos with different mass is essential to identify the best detection 
prospects for future radio telescopes like the SKA. 

Figure ^ shows the redshift-mass exclusion plot obtained by using the SKA sensitivity bound 
for SKAl LOW and SKAl MID (at 1 GHz in Band 1). For each DM halo we obtain the DM halo 
mass and redshift combination that produce the minimal SKA-detectable fluxes. DM-dominafed 
objecfs lying above the black and green curves cannot be detected with the SKAl at the given confi¬ 
dence fhreshold for (ctF) = 3 x 10^^^ cm^ s^^ Objecfs below a curve are visible fo SKAl, and the 
further below the curve they lie the greater the region of the cross-section parameter space we can 
explore through the observation of the object. For reference, the yellow dash-dotted line displays 
the curve given for {oV) = 3 x lO^^*’ cm^ s^^ and la confidence level. A few representative know 
objects (irrespective of their location in the sky) with good estimates of the DM mass are plotted in 
the Mdm — z plane for the sake of illustration of the DM search potential with the SKA. 

Dwarf galaxies, given their extreme proximity, provide an excellent test-bed for DM radio probes, 
granting access to a parameter space that extends even below the value {oV) = 3 x cm^ 

s^^. Additionally their large mass-to-light ratios and absence of strong star formation and diffuse 
non-thermal emission make them very clean sources for radio DM searches. 

Galaxies can be probed to significantly larger redshifts than the dwarf galaxies due to their larger 
DM mass, and those located in the redshift range 0.5 z ~ 1.0 provide stronger constraints. How¬ 
ever, an optimized DM search should be confined to galaxies with little background radio noise. 
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Figure 4: Exclusion plot in redshift versus halo mass based on projected SKA (at 1 GHz in Band 1) sensi¬ 
tivity data for the reference value of (crV) = 3 x 10^^^ cm ^ s^* with 30 hour integration time (black lines) 
and 1000 hour integration time (green lines). {B) = 5 /rG was adopted. Solid lines are the 1(7 sensitivity ex¬ 
clusion, dashed lines that of 2(7 and dotted lines correspond to 3(7. The yellow dash-dotted line corresponds 
to 30 hours of integration and 1(7 confidence with {oV) = 3 x 10^^** cm^ s^*. An annihilation channel bb 
is assumed with a neutralino mass of 60 GeV. Representative objects with known DM mass are shown for 
illustrative purposes of DM radio signal detection. The dSph group contains the galaxies; Draco, Sculptor, 
Fornax, Carina and Sextans. Unlabelled Galaxies are; NGC3917, NGC3949 and NGC4010. For very local 
objects the redshift is estimated from the average distance data. From Colafrancesco et al. 2014. 


making low star-formation-rate galaxies good candidates. High-z galaxies come also with the ad¬ 
vantage of observing more primitive structures with fewer sources of baryonic radio emission. 
Clusters of galaxies provide extremely good candidates in cases, such as the Bullet cluster, where 
the dark and baryonic matter are spatially separated. Our recent analysis of the ATCA observation 
of the Bullet cluster (Colafrancesco and Marchegiani 2014) indicates that deeper radio observa¬ 
tions (possible with the SKA) will be able indeed to separate the DM-induced signal from the 
CR-induced one and hence have the possibly to investigate the nature of DM particles using the 
technique here proposed. More in general, the large predicted radio fluxes due to DM annihilation 
in clusters indicate that DM-induced radio emission can be observed in radio out to large redshifts 
z « 2, again with the advantage of fewer sources of baryonic radio emission. 

2.3 Disentangling magnetic fields and Dark Matter 

Studying the magnetic properties of DM halos are crucial to disentangle the DM particle den¬ 
sity from the magnetic field energy density contributing to the expected synchrotron radio emission 
from DM annihilation. The SKA is the most promising experiment to determine the magnetic field 
structure in extragalactic sources (see Johnston-Hollitt et al. 2015), and will have the potential of 
measuring RMs toward a large number of sources allowing a detailed description of the strength, 
structure, and spatial distribution of magnetic fields in dSph galaxies, galaxies (see Beck et al. 
2015) and galaxy clusters (see Govoni et al. 2015). We stress that these measurements of the mag¬ 
netic field can be obtained by the SKA simultaneously, for the first time, with the constraints on 
DM nature from the expected radio emission. 
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3. The impact of SKA on the search for the DM nature 

Deep observations of radio emission in DM halos are not yet available, and this limits the 
capabilities of the current radio experiments to set relevant constraints on DM models. We have 
already explored a project (Regis et al. 2014a,b,c) dedicated to the WIMP search making use of 
radio interferometers, that could be considered as a pilot experiment for the next generation high- 
sensitivity and high-resolution radio telescope arrays like the SKA. For the particle DM search we 
are interested in, the use of multiple array detectors having synthesized beams of ~ arcmin size 
has a number of advantages with respect to single-dish observations. First, the large collecting area 
allows for an increase in the sensitivity over that a single-dish telescope. The best beam choice 
for the detection of a diffuse emission requires a large synthesized beam (in order to maximize the 
integrated flux), but still smaller than the source itself to be able to resolve it. A good angular reso¬ 
lution is also crucial in order to distinguish between a possible non-thermal astrophysical emission 
and the DM-induced signal, which clearly becomes very hard if the DM halo is not well-resolved. 
The possibility of simultaneously detecting small scale sources with the long-baselines of the array 
allows one to overcome the confusion limit. In the case of arcmin beams, the confusion level can 
be easily reached with observations lasting for few tens of minutes, even by current telescopes. A 
source subtraction is thus a mandatory and crucial step of the analysis. Finally, single dish tele¬ 
scopes face the additional complication related to Galactic foreground contaminations, which are 
instead subdominant for the angular scales typically probed by telescope arrays at GHz frequency. 

The limits derived from ATCA observations of 6 dSphs (Regis et al. 2014c) on the WIMP 
annihilation/decay rate as a function of the mass for different final states of annihilation/decay 
are already comparable to the best limits obtained with y-ray observations and are much more 
constraining than what obtained in the X-ray band or with previous radio observations (Spekkens 
et al. 2013, Natarayan et al. 2013). In this context, the SKA will have the possibility to explore DM 
models with cross-section values well below the DM relic abundance one (see Fig.6 in Regis et al. 
2014c). The SKAl-MID Band 1 (350 -1050 MHz) will probably be the most promising frequency 
range for the majority of WIMP models. The full SKA-2 phase will bring another factor ~ lOx 
increase in sensitivity and an extended frequency range up to at least 25 GHz. Typical values of the 
SKA sensitivity {A^ffjTsys = 2 x 10'^ m^/K) and bandwidth (300 MHz at GHz frequency) provide 
rms flux values of 30 nJy for 10 hours of integration time. This is about a 10^ factor of gain 
in sensitivity with respect to the most recent ATCA observations (Regis et al. 2014a). A further 
improvement by a factor of 2-3 can be confidently foreseen due to the larger number of accessible 
dSph satellites from the southern hemisphere. The SKA will also have the unique advantage to be 
able to determine the dSph magnetic field (via FR measurements and possibly also polarization), 
provided its strength is around the /rG level (as expected from star formation rate arguments, Regis 
et al. 2014c). This will make the predictions for the expected DM signal much more robust 
and obtainable with a single experimental configuration. The prospects of detection/constraints 
of the WIMP particle properties with the SKA will therefore progressively close in on the full 
parameter space, even in a pessimistic sensitivity case, and up to ~ TeV WIMP masses, irrespective 
of astrophysical assumptions. 

The SKA will also allow to investigate the possibility that point-sources detected in the prox¬ 
imity of the dSph optical center might be associated to the emission from a DM cuspy profile. 
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This possibility is likely only in the "loss at injection" scenario, while spatial diffusion should in 
any case flatten the e* distribution, making the source extended rather than point-like. The in¬ 
vestigation of these sources with the SKA will deserve particular attention, since we have already 
found that the WIMP scenario can fit the point-like emission with annihilation rates consistent with 
existing bounds (Regis et al. 2014c). 



Figure 5: The {oV) upper limits from 30 hour of SKA integration time for z = 0.01 at 300 MHz (top) and 
1 GHz (bottom) as the neutralino mass M'), is varied with annihilation channel bb in solid lines and Tf in 
dashed lines. A value {B) = 5 jiG was adopted. Black lines correspond to halos with mass 10*^ M©, red 
lines to 10*^ M© and green lines to 10^ M© (from Colafrancesco et al. 2014). 

The SKAl-MID Band 1 (350-1050 MHz) to Band 4 (2.8-5.18 GHz) are important to probe 
the DM-induced synchrotron spectral curvature at low-v (sensitive to DM composition) and at 
high-v (sensitive to DM particle mass), and the implementation of Band 5 (4.6-13.8 GHz) will 
bring further potential to assess the DM-induced radio spectrum. As we can see from Fig.^, the 
best frequency range to detect these radio emissions is around 1 GHz. So, the upper frequency 
regions of SKAl-MID Band 1 provides the strongest spectral candidate for probing the cross- 
section parameter space due to an optimal combination of the SKA sensitivity within this band and 
the relatively strong fluxes at these frequencies. This frequency band will also allow for an optimal 
description of the magnetic field (see Johnston-Hollitt et al. 2014). 

There are two main caveats in the forecasts for DM detection in the radio frequency band. The 
first stems from the fact that, for an extended radio emission,the confusion issue becomes stronger 
and stronger as one tries to probe fainter and fainter fluxes. Thus, the source subtraction proce¬ 
dure becomes crucial and this can affect the estimated sensitivities. The impact of this effect on 
the actual sensitivity is hardly predictable at the present time, especially for the SKA, since it will 
depend on the properties of the detected sources, the efficiency of deconvolution algorithms, and 
the accuracy of the telescope beam shape. 

The second caveat is that by bringing down the observational threshold, one can possibly start to 
probe the very low levels of possible non-thermal emission associated to the tiny rate of star forma¬ 
tion in dSph, or in galaxies and galaxy clusters. The DM contribution should be then disentangled 
from such astrophysical background. The superior angular resolution of the SKA will allow for the 
precise mapping of emissions, putatively either DM or baryonically induced, and will enable their 
correlation with the stellar or DM profiles (obtained via optical and/or kinematic measurements). 
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Early DM science can be done with a small sample of local dSphs, a small sample of nearby 
galaxies with good DM density profile reconstruction, and the Bullet cluster, observing with a 
somewhat larger beam (r:^ 1" — 10")- These objects have been already studied in radio with similar 
objectives (e.g., DM limits) and therefore provide the best science cases to prove the capabilities 
of the SKAl for the study of the nature of DM with radio observations. The implementation 
of the SKAl-MID Band 5 (4.6-13.8 GHz) will increase the ability to detect the expected high- 
frequency spectral curvature of DM-induced radio emissions, and the ability to place constraints 
on the DM cross-section and to differentiate between different annihilation channel spectra. The 
10 X increased sensitivity of SKA2-MID compared to SKAl-MID will allow us to increase the 
angular resolution by a factor 20 or to increase the sensitivity of SKA2-MID to DM-induced 
radio signals by a factor 10. The possibility to extend the frequency coverage of the SKA in its 
Phase-2 realization up to ~ 25 GHz, will allow to detect the expected high-v spectral cut-off of 
DM-induced radio emissions, and then set accurate constraints to the DM particle mass. 

4. Conclusion 

The SKA has the potential to unveil the elusive nature of Dark Matter. Its ability to resolve 
the intrinsic degeneracy (between magnetic field properties and particle distribution) of the syn¬ 
chrotron emission expected from secondary particles produced in DM annihilation (decay) will 
allow such a discovery to be unbiased and limited only by the sensitivity to the DM particle mass 
and annihilation cross-section (decay rate). The unprecedented sensitivity of the SKA to the DM 
fundamental properties will bring this instrument in a leading position for unveiling the nature of 
the dark sector of the universe. 

The information provided by the SKA can be complemented with analogous studies in other spec¬ 
tral bands, which will be able to prove the ICS signal of DM-produced secondary electrons (span¬ 
ning from /iwaves to hard X-rays and /-rays) and the distinctive presence of the //emission 
bump in the /-rays (see Fig. ||). The next decade will offer excellent multi-frequency opportunities 
in this respect with the advent of Millimetron, the largest space-borne single-dish mm. astronomy 
satellite operating in the 10^ — 10^ GHz range (optimal to prove the DM-induced SZ effect), the 
Astro-H mission operating in the hard X-rays frequency range (with the highest expected sensitiv¬ 
ity to probe the high-energy tail of the DM-induced ICS emission), and the CTA with unprecedent 
sensitivity in the energy range between a few tens GeV to hundreds TeV. 
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